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Abstract 

It has recently been proposed that violations of Lorcntz invariance or violations of the equiv- 
alence principle can be constrained from the non-observation of neutrinoless double beta decay. 
We generalize this analysis to all possible new gravitational interactions and discuss briefly the 
constraints for different cases. 



Although there is no evidence for the violation of gravitational laws, lots of work has been 
done to find out to what accuracy this is true. Many experiments have been performed to test 
the equivalence principle M for ordinary matter and to test local Lorentz invariance [Q, H]. In 
recent times there has been some effort to test these laws for the gravitational couplings of neutri- 
nos. Assuming that neutrinos of different generations have characteristic couplings to gravity with 
differing strength and that the gravitational eigenstates differ from the mass eigenstates, one can 
constrain the amount of violation of equivalence principle (VEP)in the neutrino sector from present 
neutrino oscillation results @, |5|. Similar bounds were also obtained for the amount of violation of 
local Lorentz invariance (VLI) , assuming that neutrinos of different generations have characteristic 
maximum attainable velocities || [?]]. Recently we have pointed out that in both these cases it is 
possible to constrain some otherwise unconstrained region in the parameter space from neutrinoless 
double beta decay ||. 

In this article we propose a general framework to study the effect of new gravitational interac- 
tions in the neutrinoless double beta decay. This formalism is similar to the one used in the study 
of K— system |], We classify all possible interactions as scalar, vector and tensor interactions. 
Since both the VEP and VLI are tensor interactions, it is expected that in both cases similar con- 
straints should be obtained, as observed. On the other hand, a recent string motivated violation of 
the equivalence principle a la Damour and Polyakov [11] is a scalar interaction. Thus the constraint 
in this case is of different nature than in the cases of VEP or VLI considered previously. The pos- 
sibile fifth force [0] discussed in the literature is a vector interaction and thus also has a different 
phenomenology. Our analysis can be extended to study the effects of gravitational interactions in 
neutrino oscillation experiments. 

We write down the most general lagrangian for interactions of neutrinos with scalar, vector and 
tensor fields in the weak basis [vi] following the general framework developed for the K— system 

1, 

C = GijViVj + G^Vi^Vj + G^Vi^v j)V (1) 

where i,j are generation indices and Gy, Gfj and G^J are scalar, vector and tensor fields re- 
spectively. We shall not work beyond the external field approximation. These fields have some 
restrictions coming from the symmetry properties and by discarding the total divergence expres- 
sions from the lagrangian, which have been discussed in ref. || in detail. We further assumed 
that the gravitational eigenstates could be different from the mass eigenstate as well as the weak 
eigenstate. For simplicity from now on we shall work in an two generation scenario, i, j = e,x with 
x = (J,,T, s. 
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We now can write down the Feynmann diagrams and hence the self energy matrix in the same 
way as in ref ||, from which the contribution to the effective hamiltonian can be read off in the 
weak basis, given by 

Ha = Gij + %G% Vtx + G% PliPu . (2) 

This hamiltonian is related to the effective hamiltonian in the mass and gravitational bases through 
unitary rotations 

H = U m H m U m 1 + U G H G U G 1 . (3) 

In absence of any new gravitational interactions the neutrino mass matrix in the mass basis \v\ v?\ 
is given by 

°V (4 ) 
2p 2p \ m 2 J 

and the gravitational interaction part of the hamiltonian is 

T (M G ) 2 _ l/rf 0\ , c , 

Here p denotes the momentum, / represents an unit matrix and rh the average mass, and for any 
quantity X we define SX = (Xi — X^), X = (X\ + X2)/2. a = S,V,T represents scalar, vector and 
tensor interactions respectively. 

The scalar, vector and tensor gravitational interactions can be written in the following forms 
so as to reproduce the correct dimensions of equation (||) , 

c a 2 

g { = 2a inn 

V o V 

g- = 2a v imip 
gj = 2a T iP 2 

In the absence of any gravitational interactions a a i = 0, Hq simply becomes the momentum of the 
neutrinos. Here we are interested in a single virtual neutrino propagating inside the nucleus with 
a particular momentum. For this reason we assume the momenta of both the neutrinos are p in 
the absence of any new gravitational interactions. Hence af — a\ = 5a a is a measure of the new 
gravitational interactions in the neutrino sector. To compare our result with the neutrino oscillation 
experiments we further assume, af + a^ = 0, i.e., there is no mean deviation from the gravitational 
laws and there is only a relative violation given by the measure 5a a . This approximation will reduce 
the the number of parameters so that we can compare the bounds from the neutrino oscillation 
experiments with the ones from neutrinoless double beta decay. 
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We tried to keep the above discussions as general as possible with the restrictions that we donot 
go beyond the perturbative regime. We assume that the corrections to gravity comes from interac- 
tions with some external scalar, vector or tensor fields only and there is no non-renormalizable higher 
dimensional operators which modifies gravity with inverse mass scales. Our general parametriza- 
tion has one drawback that although we are working in the gravitational basis, the masses involved 
in the expressions for gi are considered in the mass basis. This can be justified by assuming VEP 
to be a small effect. In the case of tensor interaction masses donot enter, only in the scalar and 
vector cases this problem appears. However, as we shall see, the final result for the scalar case 
comes out to be the same as the one derived from other approaches Moreover, in the case of 
some scalar interactions the gravitational basis is equal to the mass basis Q, and then this question 
will not arise. 

We shall not consider any CP violation, and hence H m and Hq are real symmetric matri- 
ces and U m and XJq are orthogonal matrices U^ 1 = U T . They can be parametrized as Ui = 

( C ° S * Sm * ), where 6i represents weak mixing angle 9 m or gravitational mixing angle 6q. 
\ — sm a i cos Oi J 

We can now write down the weak Hamiltonian H w in the weak basis, in which the charged lepton 
mass matrix is diagonal and the charged current interaction is also diagonal, as 



H 



2p 



1 ( M + M12 



\M 12 



1 



1 



pI+ — (M 1 + -M 2 M^ 



(6) 



where, 



Mi 
M 2 



UmMmU' 1 

UgMgUc 1 



and we assumed M 2 3> Mf, so that the gravitational effects are much smaller than the usual 
neutrino masses. Since no new gravitational effects have been observed so far, we use this formalism 
to constrain the parameters of the new gravitational interactions, for which this assumption is 
justified. We then obtain, 
cos29 m 



M± 



m ± 
± 



M 



2 

sin 26t 



-5m 

_ cos26» G , Sm5g a 
og m ± — - — cos 2(0g 



+ 5mg' 



cos 26 r , 



/[2(5m 2 - fh 2 )] 



12 



-5m 



5g a m 



sin 20 Q 



5mg a - 



sin 29 r , 



/2(5m 2 



m 



(7) 



1 this point will be discussed in a forthcoming article, where the dilaton-exchange gravity will be studied by the 
authors 
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to a leading order in 5g a . 

The decay rate for the neutrinoless double beta decay is given by, 

{ o^ r i = ^± Goi | MjB | 2) (8) 

where ME denotes the nuclear matrix element, Goi corresponds to the phase space factor defined 
in [|D|] and m e is the electron mass. The momentum dependence of M + must be absorbed into 
the nuclear matrix element, so that this quantity contains all the momentum dependence and the 
remaining part is estimated using zero momentum transfer approximation. Thus, if one ignores 
the nuclear matrix element, then obviously there cannot be any effect of the vector and tensor 



type gravitational interactions in neutrinoless double beta decay, which was mistaken in ref. [14]. 
As it has been discussed earlier [||, the momentum dependence of the tensor type gravitational 
interactions enters the nuclear matrix element, which then is enhanced by a factor p 2 coming 
through M + in the above expression. 

We shall now present a more detail explanation of this analysis. 

In ref [14] it is claimed that neither violations of Lorentz invariance nor violations of the equiv- 
alence principle may give sizable contributions to neutrinoless double beta decay. The argument 
discussed is the following: Taking the neutrino propagator 

e -iq{x-y)l m \ „2 

d A a — { ' a (9) 

" 2 4 221-22 V / 

™ CI ~ Qo C a + Q ' Ca 

with the standard Ov/3/3 observable (m), the neutrino four momentum q and the characteristic 
maximal velocity c a . If one would neglect now qo and m in the denominator, c a drops out and the 
decay rate is independent of c a . 

However, in || it has been shown starting from the Hamiltonian level that the propagator (or 
the Qv(3(3 observable) is changed itself violating Lorentz invariance. Since 

2 4 

H = qc a + —^ 

m (*)2 4 

with c a = I + Sv and rrS*^ 2 = vn? + 2<fc a 5v an additional contribution is obtained cx (f5v . This 
mass-like term has a q 2 enhancement and is not proportional to the small neutrino mass. This 
consideration answers also the frequently asked question "What is the source of lepton number 
violation?" in this mechanism. Comparable to a usual mass term, which can be both of Majoran 
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type as well as Dirac type the mass-like term 2(fc a 5v can be of Majorana type and act as the 
source of lepton number violation in this context. 

We shall now discuss the three different cases of scalar, vector and tensor interactions and their 
phenomenology. In the case of tensor interaction the constraint has already been discussed in ref. 
||]. The violation of local Lorentz invariance and the violation of the equivalence principle both 
fall under this category (their equivalence has been pointed out elsewhere [§]). In both these cases 
the effect of the new gravitational interactions have quadratic momentum dependence. 

In case of the tensorial gravitational interaction we have, g T = 2a T p 2 = and 5g T = 25a T p 2 . 
In particular, for the violation of the equivalence principle we substitute 5a T = 4:5g4> (following the 
notation of ref pi), where <fi is the Newtonian gravitational potential on the surface of earth. On 
the other hand, for the violation of the local Lorentz invariance, we substitute instead, 5a T = 2Sv. 
In both these cases bounds were given in ref ||. 

To give a bound on tensorial gravitational interactions in the small mixing region (including 
&v ~ G m ~ 0) conservatively (m) ~ was assumed. It was also assumed that 5m < m, and 
thus ||j may be neglected. Due to the p 2 enhancement the nuclear matrix elements of the mass 
mechanism have to be replaced by ^ • (M' F — M' GT ) with the nuclear radius R and the proton 
mass m p , which have been calculated in |j~5[j . Inserting the recent half life limit obtained from the 
Heidelberg-Moscow experiment [16|, a bound on the amount of tensorial gravitational interactions 
as a function of the average neutrino mass fh was given H. It should be stressed also that the 



GENIUS proposal of the Heidelberg group [17] could improve these bounds by about 1-2 orders of 
magnitude. 

For the vector type gravitational interactions there is a linear momentum dependence. In this 
case, g v = 2a v pm = and 5g v = 28a v pm. The fifth force, as discussed by Fishbach et al |l^] in 
the context of K— physics is a vector type gravitational interaction. Since no studies of this type of 
forces exist for neutrino oscillation experiments, with which neutrinoless double beta decay results 
could be compared, we shall not study this case. 

A similar generic structure was considered in a recent analysis of the atmospheric neutrino 
anomaly [JOI], where they used the power of momentum dependence as a parameter. From their 
analysis it becomes apparent that the atmospheric neutrino anomaly may not be explained by 
either tensorial or vectorial gravitational analysis alone [O]. 

Recently it has been argued by Damour and Polyakov [|ll]] that string theory may lead to a new 
scalar type gravitational interaction through interaction of the dilaton field and subsequently its 
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consequence to neutrino oscillation has been studied Damour and Polyakov have shown that 
the massless dilaton interaction modifies the gravitational potential energy and there is an addi- 
tional contribution from an spin-0 exchange, which results in a scalar type gravitational interaction 
| |1 lfl . The resulting theory is of scalar-tensor type with the two particle static gravitational energy 

V(r) = -G N m A m B (l + ctA<XB)/r, (11) 

where Gat is Newton's gravitational constant and ay denotes the couplings of the dilaton field <j) to 
the matter field ipj, leading to a gravitational energy of 

L = mjajipjifjj(ft. (12) 



Thus the modified effective mass matrix of the neutrinos are now given by [19] 

my 1 = m — ma<j) c (13) 

where, the classical value of the dilaton field <p c = 4>N^ext is characterised by the a value of the bulk 
matter producing it and for a static matter distribution proportional to the Newtonian potential 

<t>N- 

The effective mass squared difference 

Am ( *» 2 = -2m 2 <j) N a ext 5a (14) 

(for almost degenerate masses m) gives rise to neutrino oscillation. The corresponding effect for 
Qvf3f3 decay is obtained by replacing 5g s = 25a s vn? (for almost degenerate mass m\ ~ ?ti,2 ~ m). 
Comparing the arguments in the oscillations propabilities we get 

M+=m + ma ex t^NSa COS ^ 0G \ (15) 

In this case, it is difficult to obtain any bound from neutrino experiments since for a ex t only 
upper bounds exist. However, to get an idea of the constraints which can come from neutrino 



experiments in the future if a ex t is known, according to ref. 1£] we assume (pN = 3 • 10 -5 , 
a ex t = VlO -3 and m = 2.5 eV (as an upper bound obtained from tritium beta decay experiments 
p0|| ). In this case the quantity 5a is not constrained from neutrinoless double beta decay. 

In summary, we presented a general formalism for the study of effects of new gravitational 
interactions in neutrinoless double beta decay, which allows to constrain the amount of violation of 
the gravitational laws. Various scenarios discussed in the literature have been analyzed as special 
cases of the present formalism. 
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